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Abstract 



Q ■ The antiferromagnetic ground state of the half-filled Hubbard model with the doubly 

cn . degenerate orbital has been studied by using the slave-boson mean-field theory which 

> 

C<S ' was previously proposed by the present author. Numerical calculations for the simple 

O' 

(N 

Q ■ cubic model have shown that the metal-insulator transition does not take place except at 

the vanishing interaction point, in strong contrast with its paramagnetic solution. The 
energy gap in the density of states of the antiferromagnetic insulator is much reduced by 

C . the effect of electron correlation. The exchange interaction J plays an important role in 

O 

the antiferromagnetism: although for J = the sublattice magnetic moment m in our 

theory is fairly smaller than tt^hfa obtained in the Hartree-Fock approximation, m for 
J/U > 0.2 {U: the Coulomb interaction) is increased to become comparable to mnpA- 
Surprisingly, the antiferromagnetic state is easily destroyed if a small, negative exchange 
interaction {J/U < —0.05) is introduced. 

PACS No. 71.27+a, 71.30.+h, 75.10.Lp 



I. INTRODUCTION 

Much progress has been made in our theoretical understanding on the effect of elec- 
tron correlation in systems such as transition metals and high-Tc materials. Most of the 
theoretical studies have been made for the single-band Hubbard model (SHM)^~^ for its 
simplicity. Actual systems, however, inevitably have the orbital degeneracy. It is neces- 
sary to investigate the role of the orbital degeneracy and the effect of Hund-rule coupling 
due to the exchange interaction for a better understanding on strongly correlated systems. 

In the last few years the Hubbard model with orbital degeneracy has been exten- 
sively studied by using various methods such as the Gutzwiller approximation (GA),"^~^ 
the slave-boson theory^'®, the dynamical mean-field approximation,^ and the projective 
self-consistent method.^'' The original GA proposed by Gutzwiller and Chao^^ was refor- 
mulated in Refs.4-6. In a previous paper'^ (referred to as I hereafter), the present author 
developed the slave-boson functional-integral method for the Hubbard model with an arbi- 
trary, orbital degeneracy, by employing the method proposed by Dorin and Schlottman^^ 
for the Anderson lattice model. Fresard and Kotliar^ developed an alternative slave-boson 
functional integral method. These slave-boson theories are the simple generalization of 
the Kotliar-Ruckenstein theory for the SHM^^ to that for the degenerate Hubbard model, 
and their saddle-point approximation is equivalent to the GA.'^^^'^^ In I we have studied 
the metal-insulator (MI) transition of the doubly degenerate Hubbard model (DHM) in 
the paramagnetic state. The MI transition takes place when the interaction strength is 
increased, ^'^'^'^'^ just as in the case of the SHM.^^ This MI transition is shown to become 
the first-order one in the half-filled case when the exchange interaction is included. ^'^ 

We should, however, remind the following facts having been established for the half- 



filled SHM: 

(a) In the paramagnetic (P) state, the MI transition is realized for d = 1,2,3 and oo in 
the GA, but not for rf = 1, 2 and 3 in the advanced theory going beyond the GA.^^^^^ 

(b) In the antiferromagnetic (AF) state, the MI transition occurs in neither d = 1, 2, nor 
3^^ in the advanced theory, nor in d = oo even within the GA.^^~^^ 

These facts suggest that it is indispensable to take into account the antiferromagnetic 
state in discussing the MI transition in DHM. 

One of the advantages of the slave-boson functional integral method over the GA 
is that it has the wider applicability than the GA. For example, we can deal with the 
system with the complicated magnetic structures such as the antiferromagnetic state, by 
using the Green's function formalism. We will study in this paper, the antiferromagnetic 
state of the DHM by employing our slave-boson mean-field theory,^ in order to clarify the 
above-mentioned issue relevant to the MI transition and the roles of the degeneracy and 
the exchange interaction. 

The paper is organized as follows: In the next Sec. II, we present a basic formulation 
of our slave-boson saddle-point approximation to deal with the antiferromagnetic state in 
the DHM, after briefly reviewing I. Numerical calculations for the simple-cubic lattice are 
presented in Sec. HI. Section IV is devoted to conclusion and supplementary discussion. 



II. FORMULATION 

We adopt the Hubbard model with the arbitrary, orbital degeneracy D, whose Hamil- 



tonian is given by 

cr ij mm' i {rn,(T)j^{m' ,a') 

where Qmo- is an annihilation operator of an electron with an orbital index m and spin 
o"(=|', I) on the lattice site i. The electron hopping is assumed to be allowed only between 
the same sub-band: t^"^ = tjj (5mm', for a simplicity. The on-site interaction, [/^'^,, is given 
by 

U^m' = Uo = U ioTm = m',(Tj^a', (2) 

= u^ = u-2j ioTm^m',a^ a', (3) 

= U2 = U-3J foTm^m',a = a', (4) 

where U and J are Coulomb and exchange interactions, respectively. 

In 1 we employed the boson opertor intoduced by Dorin and Schlottman,^^ and used 
the static approximation to get the functional integral representation of the partition 
function given by^ 

Z = f D^ I Du I Dm I Dn I Il'j^^Db'-^^ exp {-(3^) , (5) 



with 



e ^* = exp —(3 



7 y / \ >,irn''''irn '^im'^im) ' ^C 



Tr exp{-pHeff), (6) 



2D 

/ y / y / y / y ^ m,m,' "i;mi CTj , . . . . ,»7i^(T^ "j;miCTi,....,m^(T£- \' ) 

i 1=2 miai,....,miai (mcr,m'a') 

The effective Hamiltonian is given by 



a ij m (Tim 



with 






2D 

\J(^i Pima + VPi^ ^f^ + H (^^ 

£=3 



(^-1) . j^W 



(9) 



In Eqs.(5)-(9), D^ = HimCim et al, ^im (i^im) is the exchange (charge) field in the subband 
m at the site i, and mj^ (^im) is the magnetic moment (electron number). The forth 
summation in eq.(7) is performed over a pair of indices {ma,m'a') with [ma) ^ (mV) in 
the configuration: {miO"i,m20'2j ■■■,'miai} occipied by £ electrons. The introduced boson 



operator, b\ : 



^i %\m\a\,....,'miiuii 



(10) 



projects to the configuration of I electrons with pairs of orbital and spin indices {mo"}. 
Its full contraction, (6^^ ■ 6^ ), and partial contraction, (6^- ■o'l^ma-, are defined by 






E ^1 



(^) 
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rri\(7\^....^miGi ^z;rn,io"i,....,m^cr^ "> 



:iil 



m\G\^....^miGi 



\^l ^l jrrinUn — / , 

mio-i,m20-2,....,m„-io-„_i,m„+io-m+ilv-.m-£0-£ 






(12) 



The empty state (cj) and singly occupied one with a cr-spin electron (pimo-) are expressed 
in terms of n,„, , m,„, and 6,- for £ > 2 as''' 



nm; ""im "'iJ-'-i "j 



,(0) ,(0)> 



2D 



e, = {hT ■ b^ = 1 - E ^- + E E [(^ - 1)/^] (^^ ^ ■ b. 

m mcr £=2 



j Jmai 



2D 



Pima = (b- ■ hi ')ma = (Uim + (Tmim)/2 - ^{hl ■ b, 



i ) ma ■ 



(13) 



(14) 



(.=2 



The expression for the functional integral given by eqs.(5)-(14) is a generalization of the 
single-band model to the degenerated-band model, ^^'' and it has a transparent physical 



meanmg. 
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When we apply our slave-boson functional integral method developed in I to the DHM 
{D = 2), the functional integral becomes 

Z = JDi JDm I Dv JDn f Ddo j Ddi fu^Dd^ fll^DU fof ea:p[-/3(<l'o + $i + $2)], 

(15) 
$0 = E Wd>x^ + 2Uid,i + f/2(rf.T + dii) + 2(f/o + f/i + U2)m + U^ + fi)] , (16) 

i 

with 

$1 = E [iirnmim + (/^ - Vim)nirn] , (17) 

^^= I de f\e) (-I/tt) Jm Tr In G{e). (18) 

In Eqs.(15)-(18), f{e) is the Fermi-distribution function, and rfj, tj and /j denote the states 
with double, triple and quadruple occupations, respectively. In particular for doubly 
occupied states, we take into account the three kinds of configurations: rfjo for a pair of 
electrons on the same orbital with opposite spin, dn on the different orbital with opposite 
spin, and rfjo- on the different orbital with same spin a. 

The one-particle Green function, G'(e), in Eq.(18) is expressed by 

G{e) = {s-H,ff)-\ (19) 

where the effective Hamiltonian, -ffe//, is given by 

(T ij m aim 

the band-narrowing factor, g^^, being given by 

9ma ^ima ^jmai V^-*-/ 



with 



.(22) 



^'^ {riim + crniimy/'^ {2 - riim - o-rriimy/'^ 

ei = l- 2nim + 2dio + 2dii + di^ + diy + 4(ti| + Uy) + 3/^, (23) 

Pima = {riim + crrriira) /2 - {dio + dii + di„) - 2 tj^ - t^,^ - /j. (24) 

In order to discuss the antiferromagnetic (AF) state, we divide the crystal into two 
sublattices, A and B. We assume that for the AF wave vector Q, the relation: Ek+q = —£k 
holds where Ek is the Fourier transform of the transfer integral, tij. We take $,im in Eq.(20) 
as the staggered field given by $,im = C (—0 ioi i E A {i E B), the exchange fields in the 
two subbands being assumed to be the same. The magnitude of ^ will be determined by 
the variational condition, as will be shown shortly (Eq.(31)). 

Since the effective transfer integral in Eq.(20) is expressed as a product form: 
ZimatijZjma, we cau express the one-electron Green function in terms of the locators defined 

by23) 

where rima = (-^imo-)^ = ^Ama and rsma for i G A and i E B, respectively. After a simple 
calculation, we get $2 given by 

^2 = Jde fie) (I/tt) Im Y. ^^ (?L [XAm.{e)XBm.{e) - £^]) , (26) 

mkcF 

where the band-narrowing factor, qma-, is given by 

DeCaUSe Zjs^xaa Zj^fYi^fj ZfYifj anQ Tj^-u^j ^Bm—a ^mu' 



The mean-field free energy is obtained from the saddle-point values of the integration 
variables for which the variational conditions yield the following simultaneous equations: 

(T 
(T 

fJ'-l^im + Yl Rima{drima/dnim) = 0, (30) 

a 

2f/o + Y. Rrma{dr,^Jdd,^) = 0, (32) 

mcr 

2Ui + J2 R^mA^r■^m./^da) = 0, (33) 

ma 

U2 + J2 Rima' (drima'/ddi^) = 0, (34) 

mcr' 

2(f/o + f/i + f/2) + E ^^--' (-^rim^'Z-^t..) = 0, (35) 

m.a' 
2{Uo + f/l + f/2) + E R^mA^nma/^f^) = 0, (36) 

In Eqs.(28)-(36), Rima and n^m^ are given by 

Rima = d^2/drima = j ^S f {s) (-I/tt) Im [{VLrna/rima) Fo(fim^)] , (37) 

Uima = de f{e) Pima{e), (38) 



where the local densities of states at the site belonging to A and B sublattices are expressed 

by 

Pimaie) = (-I/tt) Jm [KAma{£)/rAma] (^ e A), 

= (-i/tt) Im [KBrna{e)/rBm.] (^ G B), (39) 



with 



KAna{e) = [XB^.{e)lXA^a{e)Y^'F^{n 



1/2 



KBnai^) — [-^Amcr(^)/-^Bmo-(^)] -fo(^mcr), 



1/2 



F,{e) = jde'p,{e')/{e-e') 
p^{e) being the unperturbed density of states.^'' 



1/2 



(40) 
(41) 
(42) 
(43) 



III. NUMERICAL CALCULATIONS 

Numerical calculations have been performed for the simple-cubic model with nearest- 
neighbor hoppings t. Input parameters for our calculations are the non-interacting den- 
sity of states, Pq{s), the Coulomb and exchange interactions, U and J, and the number 
of electrons per sub-band, n, which is unity for the half-filled case. We employed the 
approximate, analytic expression for po(^) of the simple-cubic lattice, given by^^ 



Po{e) = A 
= A 
= 



9-iu' 



9-uj' 



1/2 



1/2 



c 



B 



l-u' 



1/2 



for I cij |< 1, 



l-{\u\-2f 



1/2 



for 1 <| u; |< 3, 



for I a; I > 3, 



(44) 



where a; = e/2t, A/2t = 0.101081, B/2t = 0.128067 and C/2t = 0.02. The energy 
and the interactions are hereafter measured in units of a half of the total band width, 
W/2 = 6 t = 1. The ground-state energy without interactions calculated by using Eq.(44) 
is Eq = —0.3349, which is in good agreement with the exact value of —0.3341.^^ Since the 
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relations: e = / and p^r = t^ hold for the half-filled case, we have to self-consistently solve 
Eqs.(28)-(43) for nine quantities: m,C,,do,di,d^,di,ti^,ti and /, by using the Newton- 
Rapson method. We performed the integrations given by Eqs.(37) and (38) with the use of 
the contour integral along the complex energy axis,^^ in order to reduce the computational 
time. 

A. J = Case 

We firstly show the calculated results for the vanishing exchange interaction (J = 
0), for which do and di are equivalent. Figure 1 shows the sublattice magnetization, 
m, as a function of U. The antiferromagnetic state is realized for an infinitesimally 
small interaction. The magnetic moment increases with increasing U and asymptotically 
approaches the saturated value of 1.0 /i^ as [/ ^ oo. Because of large fluctuations, m 
in GA is much reduced by more than 50 % than that in the Hartree-Fock approximation 
(HFA) at f/ < 1. 

Figure 2 shows the spin- dependent local densities of states for U = 1.0. They have 
traces of the van Hove singularity of the simple-cubic density of states and clear en- 
ergy gaps characteristic of the antiferromagnetic insulator. The energy gap in the GA 
(Aga = 0.114) is much reduced compared with that in the HFA (Ahfa = 0.764). The 
U dependence of the energy gap is plotted in Fig.l. Both Aga and Ahfa increase with 
increasing U because A = 2 ^, ^ being the staggered exchange field. The ratio defined by 
a = Aga/Ahfa = ^ga/Chfa is unity in the limits of f/ ^ and f/ -^ oo, and it has a 
broad minimum of a = 0.12 at f/ ~ 0.6, above which a again increases: a = 0.15, 0.29 
and 0.54 for U = 1.0, 1.5 and 2.0, respectively. 
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The band-narrowing factor, which becomes gmt — Imi = q in our half-filled model, 
is shown as a function of U in Fig.3. In the P state, q monotonously decreases with 
increasing U as shown by the dotted curve, and it vanishes at f/ = f/c = 12^0 = 4.019 
where Uc denotes the critical interaction for the MI transition. ^'^'^ On the contrary, the 
U dependence of q in the AF state is quite different from that in the P state. When U is 
increased from the zero value, q of the AF state gradually departs from that of the P state, 
and it has the minimum value of 0.837 at [/ = 1.4, above which q increases again. The 
effect of electron correlation on the band-narrowing factor is not considerable although 
its effect on the energy gap (or the exchange field) is significant. 

The [/-dependence of the occupancies is shown in Fig. 4. At U = all the occupancies 
are 0.0625 (= 1/2^). When U value is increased, only d^ considerably increases, approach- 
ing unity for U = oo: do (= di) and t| have small peaks at t/ ~ 1 but decrease for larger 
U. 

The U dependence of the ground-state energies, E, is shown in Fig. 5. The ground- 
state energy of the AF state (Eaf) calculated by the GA is not only lower than that of 
the P state (-Ep) obtained by the GA but also lower than that of the AF state calculated 
by the HFA. The difference: AE = ii^AF(GA) — EafIUFA) expresses the energy gain by 
including the effect of fiuctuations, and its maximum value is —0.056 at f/ = 0.95. The 
HFA for the Neel state is a good description of the half-filled DHM in the limit of f/ = oo. 



B. Finite J Case 



Next we introduce the exchange interaction, J, into our calculation. Figure 6 shows the 
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sublattice magnetization as a function oi U + J for various choices of the ratio: J/U = 0.1, 
0.2 and 0.3. Note that the magnetization in the HFA is universal when it is plotted against 
U + J because its exchange field is given by .^hfa = (1/2) (f^ + <^)^hfa- As the value of 
J/U is increased, the sublattice magnetization is increased as expected. This fact is more 
clearly seen in Fig. 7, where the sublattice magnetization and the band-narrowing factor 
for U = 1.0 are plotted as a function of J. The sublattice magnetization of the GA, 
particularly near J = 0, is much increased when the J value is increased, although such 
an increase in m is realized also in the HFA result, but very small. 

The interaction dependence of the band narrowing-factor is shown in Fig.8. It was 
recently pointed out^'^ that, when J is finite, the first-order MI transition is realized in 
the P state, as is shown by dotted curves; it occurs at f/ + J = 2.21, 1.95 and 1.83 for 
J/U = 0.1, 0.2 and 0.3, respectively. Our calculation shows that the situation is quite 
different in AF state: q decreases only slightly and never vanishes. The minimum values 
of q are 0.939, 0.965 and 0.975 for J/U = 0.1, 0.2 and 0.3, respectively. 

Figure 9 shows the occupancies as a function of f/ + J in the typical case of J/U = 0.1. 
For J > the degeneracy between do and di is removed and we get do < di. When the 
interaction is increased, only d^ has an appreciable value at U + J > 0.5, as in the case 
shown in Fig. 4. 

Figure 10 shows the ground-state energies, Eaf{GA), Ej<^p(B.FA) and Ep{GA), as a 
function of the interaction. We realized that -Eaf(GA) is the lowest among the three for 
any U investigated. The maximum difference of AE is —0.029 at f/ = 0.7. 

IV. CONCLUSION AND DISCUSSION 
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To summarize, we have studied the antiferromagnetic ground state in the DHM, em- 
ploying our slave-boson mean- field theory/ Numerical calculations have shown that the 
MI transition does not take place in the antiferromagnetic solution for the half-filled DHM 
except at the vanishing interaction point, which arises form a peculiarity due to the perfect 
nesting in the model. This is in contrast with the result in its paramagnetic solution, ^'^ 
but is the same as the half-filled SHM,^^ as was discussed in the Introduction. Except at 
U = J = 0, the stable state is the antiferromagnetic insulator, whose energy gap is much 
reduced by electron correlation. 

It is worth to make a brief comparison between the results of the DHM and SHM. 
Dashed lines in Figs.l, 3, 6 and 8 show the interaction dependence of m, A and q of 
the SHM.^^ When we compare these results with the corresponding ones of the DHM, 
we notice that both the results are very similar provided the exchange interaction is not 
small; J/U > 0.2. When J is small, however, m, A and q in the DHM are fairly smaller 
than those in the SHM. Figure 7 shows that the exchange interaction effectively works to 
increase the magnitude of sublattice moment in the DHM. 

In order to more investigate the role of the exchange interaction on antiferromagnetism 
in the DHM, we have repeated a numerical calculation for the negative J, although J is 
conventionally taken to be positive. We notice in Fig. 7 that when the negative exchange 
interaction is included, the sublattice magnetization for ^ = 1.0 is considerably reduced 
and it disappears for J/U < —0.05. Figure 11 shows the sublattice magnetization and 
the band narrowing factor as a function oi U + J for J/U = —0.02. When the interaction 
is increased, the magnetization first increases at U + J < 1.2 but decreases at larger 
interaction. Surprisingly the antiferromagnetic state disappears at f/ + J > 2.1, which is 
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in strong contrast with the HFA result shown by the dotted curve. Figure 12 shows the 
interaction dependence of the occupancies. The behavior of the double occupancies in the 
negative J case is rather different from that in the positive J case shown in Fig. 9. When 
J is negative, the doubly occupied state with the opposite spin between the different 
subbands is less favorable than that within the same subband. Furthermore the triplet 
state expressed by d^ (or rfj) becomes less stable than the singlet state expressed by d^ or 
(ii, which works to suppress the antiferromagnetism. We get d^ > di > d^ {= d^) m. the 
paramagnetic state at f/ + J > 2.1. The first-order MI transition occurs at f/ + J = 2.68. 
A material with the anti-Hund-rule coupling (J < 0) would show the unusual behavior if 
it exists. 

The model Hamiltonian adopted in our study (Eq.(l)) is relevant to systems with 
partially filled narrow degenerate bands. A typical example is V2O3, which is an antifer- 
romagnetic insulator (AFI) in the ground state and which shows the MI transition between 
AFI, paramagnetic metal and paramagnetic insulator as a function of the temperature, 
the pressure and/or the chemical substitution. This phase diagram can be qualitatively 
understood with the SHM.^^'^^ The degenerate Hubbard model has much variety than 
the SHM because it has an additional, orbital degree of freedom. It may show the orbital 
ordering besides the spin ordering. The stability of the different phases in the degen- 
erate band depends not only on the values of the various interactions included in the 
model Hamiltonian but also on the temperature. It would be interesting to investigate 
the temperature-interaction phase diagram of the DHM, by generalizing our approach^^ in 
which the effects of electron correlation and thermal spin fluctuations are properly taken 
into account. 
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Figure Captions 

Fig. 1 The sublattice magnetization, m, and the energy gap, A, as a function of U of 
the DHM with J = in the GA (sohd curve) and in the HFA (dotted curve), the results 
of the SHM being shown by the dashed curve (Ref.21). 

Fig. 2 The spin-dependent local densities of states for U = 1.0 and J = 0.0: |-spin (solid 
curve) and |-spin (dashed curve) components in the GA, and |-spin (dot-dashed curve) 
and |-spin (dottedcurve) components in the HFA. 

Fig. 3 The band-narrowing factor, q, as a function of U of the DHM with J = in the 
AF state (solid curve) and in the P state (dotted curve), the result of the SHM in the AF 
state being shown by the dashed curve (Ref.21). 

Fig. 4 The occupancies as a function of U with J/U = 0, the result of d^ divided by a 
factor of ten being plotted by the dot-dashed curve. 

Fig. 5 The [/-dependence of the ground-state energies E of the AF state with J = 
in the GA (solid curve) and in the HFA (dashed curve), and their difference: AE = 
Eaf{GA) — Exy{HFA). The result of the P state in the GA is shown by the dotted 
curve. 

Fig. 6 The sublattice magnetization, m, as a function oiU + J oi the DHM with various 
J/U values in the GA (solid curve) and in the HFA (dotted curve), the result of the SHM 
being shown by the dashed curve (Ref.21). 
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Fig. 7 The sublattice magnetization for U = 1.0 as a function of J in the GA (sohd 
curve) and in the HF (dotted curve), the result of the band narrowing factor being also 
plotted (chain curve). 

Fig. 8 The band-narrowing factor, q, as a function of f/ + J of the DHM with various 
J/U values in the AF state (solid curve) and in the P state (dotted curve), the result of 
the SHM in the AF state being shown by the dashed curve (Ref.21). 

Fig. 9 The occupancies for J/U = 0.1 as a function oi U + J, the result of d^ divided by 
a factor of ten being plotted by the chain curve. 

Fig. 10 The interaction dependence of the ground-state energies of the AF state with 
J/U = 0.1 in the GA (solid curve) and in the HF (dashed curve), and their difference: 
AE = Eaf{GA) - Eaf{HFA). The result of the P state in the GA is shown by the 
dotted curve. 

Fig. 11 The sublattice magnetization, m, as a function oi U + J with J/U = —0.02 in 
the GA (solid curve) and in the HFA (dotted curve), the band narrowing factor, q, being 
shown by the chain curve. 

Fig. 12 The occupancies for J/U = —0.02 as a function of f/ + J. 



